Dynamic Global Vegetation Model Research Articles

Climate change critically affects the status of the land-system change planetary boundary

The planetary boundaries framework defines a safe operating space for humanity. To date, these boundaries have mostly been investigated separately, and it is unclear whether breaching one boundary can lead to the transgression of another. By employing a dynamic global vegetation model, we systematically simulate the strength and direction of the effects of different transgression levels of the climate change boundary (using climate output from ten phase 6 of the Coupled Model Intercomparison Project models for CO2 levels ranging from 350 ppm to 1000 ppm). We focus on climate change-induced shifts of Earth’s major forest biomes, the control variable for the land-system change boundary, both by the end of this century and, to account for the long-term legacy effect, by the end of the millennium. Our simulations show that while staying within the 350 ppm climate change boundary co-stabilizes the land-system change boundary, breaching it (>450 ppm) leads to critical transgression of the latter, with greater severity the higher the ppm level rises and the more time passes. Specifically, this involves a poleward treeline shift, boreal forest dieback (nearly completely within its current area under extreme climate scenarios), competitive expansion of temperate forest into today’s boreal zone, and a slight tropical forest extension. These interacting changes also affect other planetary boundaries (freshwater change and biosphere integrity) and provide feedback to the climate change boundary itself. Our quantitative process-based study highlights the need for interactions to be studied for a systemic operationalization of the planetary boundaries framework.

Conflicting Changes of Vegetation Greenness Interannual Variability on Half of the Global Vegetated Surface

AbstractChanges in the interannual variability (IAV) of vegetation greenness and carbon sequestration are key indicators of the stability and climate sensitivities of terrestrial ecosystems. Recent studies have examined the changes in the vegetation IAV using atmospheric CO2 observations and dynamic global vegetation models (DGVMs), however, reported different and even contradictory IAV trends. Here, we investigate the changes in the IAV of vegetation greenness, quantified as coefficient of variability (CV), over the past few decades based on multiple satellite remote sensing products and DGVMs. Our results suggested that, on half of the global vegetated surface (mostly in the tropics), the CV trends detected by different satellite remote sensing products are conflicting. We found that 22.20% and 28.20% of the global vegetated surface (mostly in the non‐tropical land surface) show significant positive and negative CV trends (p ≤ 0.1), respectively. Regions with higher air temperature and greater aridity tend to have increasing CV trends, whereas greater vegetation greening trend and higher nitrogen deposition lead to smaller CV trends. DGVMs generally cannot capture the CV trends obtained from satellite remote sensing products, while the inconsistency among satellite remote sensing products is likely caused by their process algorithms rather than the sensors utilized. Our study closely examines the changes in the IAV of global vegetation greenness, and highlights substantial uncertainty when using satellite remote sensing to study the response of terrestrial ecosystems to climate change.

Biospheremetrics v1.0.2: an R package to calculate two complementary terrestrial biosphere integrity indicators – human colonization of the biosphere (BioCol) and risk of ecosystem destabilization (EcoRisk)

Abstract. Ecosystems are under multiple stressors, and impacts can be measured with multiple variables. Humans have altered mass and energy flows of basically all ecosystems on Earth towards dangerous levels. However, integrating the data and synthesizing conclusions is becoming more and more complicated. Here we present an automated and easy-to-apply R package to assess terrestrial biosphere integrity that combines two complementary metrics. (i) The BioCol metric that quantifies the human colonization pressure exerted on the biosphere through alteration and extraction (appropriation) of net primary productivity.(ii) The EcoRisk metric that quantifies biogeochemical and vegetation structural changes as a proxy for the risk of ecosystem destabilization. Applied to simulations with the dynamic global vegetation model LPJmL5 for 1500–2016, we find that large regions presently (period 2007–2016) show modification and extraction of >20 % of the preindustrial potential net primary production. The modification (degradation) of net primary production (NPP) as a result of land use change and extraction in terms of biomass removal (e.g., from harvest) leads to drastic alterations in key ecosystem properties, which suggests a high risk of ecosystem destabilization. As a consequence of these dynamics, EcoRisk shows particularly high values in regions with intense land use and deforestation and in regions prone to impacts of climate change, such as the Arctic and boreal zone. The metrics presented here enable spatially explicit global-scale evaluation of historical and future states of the biosphere and are designed for use by the wider scientific community, being applicable not only to assessing biosphere integrity but also to benchmarking model performance. The package will be maintained on GitHub and through that we encourage its future application to other models and data sets.

Inclusion of bedrock vadose zone in dynamic global vegetation models is key for simulating vegetation structure and function

Abstract. Across many upland environments, soils are thin and plant roots extend into fractured and weathered bedrock where moisture and nutrients can be obtained. Root water extraction from unsaturated weathered bedrock is widespread and, in many environments, can explain gradients in vegetation community composition, transpiration, and plant sensitivity to climate. Despite increasing recognition of its importance, the “rock moisture” reservoir is rarely incorporated into vegetation and Earth system models. Here, we address this weakness in a widely used dynamic global vegetation model (DGVM; LPJ-GUESS). First, we use a water flux-tracking deficit approach to more accurately parameterize plant-accessible water storage capacity across the contiguous United States, which critically includes the water in bedrock below depths typically prescribed by soil databases. Secondly, we exploit field-based knowledge of contrasting plant-available water storage capacity in weathered bedrock across two bedrock types in the Northern California Coast Ranges as a detailed case study. For the case study in Northern California, climate and soil water storage capacity are similar at the two study areas, but the site with thick weathered bedrock and ample rock moisture supports a temperate mixed broadleaf–needleleaf evergreen forest, whereas the site with thin weathered bedrock and limited rock moisture supports an oak savanna. The distinct biomes, seasonality and magnitude of transpiration and primary productivity, and baseflow magnitudes only emerge from the DGVM when a new and simple subsurface storage structure and hydrology scheme is parameterized with storage capacities extending beyond the soil into the bedrock. Across the contiguous United States, the updated hydrology and subsurface storage improve annual evapotranspiration estimates as compared to satellite-derived products, particularly in seasonally dry regions. Specifically, the updated hydrology and subsurface storage allow for enhanced evapotranspiration through the dry season that better matches actual evapotranspiration patterns. While we made changes to both the subsurface water storage capacity and the hydrology, the most important impacts on model performance derive from changes to the subsurface water storage capacity. Our findings highlight the importance of rock moisture in explaining and predicting vegetation structure and function, particularly in seasonally dry climates. These findings motivate efforts to better incorporate the rock moisture reservoir into vegetation, climate, and landscape evolution models.

A new temperature–photoperiod coupled phenology module in LPJ-GUESS model v4.1: optimizing estimation of terrestrial carbon and water processes

Abstract. Vegetation phenological shifts impact the terrestrial carbon and water cycle and affect the local climate system through biophysical and biochemical processes. Dynamic global vegetation models (DGVMs), serving as pivotal simulation tools for investigating climate impacts on terrestrial ecosystem processes, incorporate representations of vegetation phenological processes. Nevertheless, it is still a challenge to achieve an accurate simulation of vegetation phenology in the DGVMs. Here, we developed and implemented spring and autumn phenology algorithms into one of the DGVMs, LPJ-GUESS. The new phenology modules are driven by temperature and photoperiod and are parameterized for deciduous trees and shrubs by using remotely sensed phenological observations and the reanalysis data from ERA5. The results show that the LPJ-GUESS with the new phenology modules substantially improved the accuracy in capturing the start and end dates of growing seasons. For the start of the growing season, the simulated RMSE for deciduous trees and shrubs decreased by 8.04 and 17.34 d, respectively. For the autumn phenology, the simulated RMSE for deciduous trees and shrubs decreased by 22.61 and 17.60 d, respectively. Interestingly, we have also found that differences in the simulated start and end of the growing season also alter the simulated ecological niches and competitive relationships among different plant functional types (PFTs) and subsequentially influence the terrestrial carbon and water cycles. Hence, our study highlights the importance of accurate phenology estimation to reduce the uncertainties in plant distribution and terrestrial carbon and water cycling.

A comprehensive evaluation of hydrological processes in a second‐generation dynamic vegetation model

AbstractThe global water and carbon cycles are greatly influenced by terrestrial vegetation, making trustworthy representations of dynamic biosphere–hydrosphere interactions a crucial component of both ecosystem and climate models. This paper comprehensively evaluates the hydrological performance of a leading dynamic global vegetation model Lund‐Potsdam‐Jena General Ecosystem Simulator (LPJ‐GUESS), using a broad range of the latest available global observation‐based gridded datasets that cover the main components of the hydrological cycle. Overall, we find that the hydrological components modelled by LPJ‐GUESS agree well with global gridded datasets of runoff, evapotranspiration and surface soil moisture, though there are discrepancies in some regions and periods. Furthermore, LPJ‐GUESS accurately captures both inter‐ and intra‐annual variations of runoff in most regions and catchment areas, including the Danube, Murray, Yangtze, Yenisei and Nile basins. Total evapotranspiration modelled by LPJ‐GUESS agrees closely with the evapotranspiration estimates of the Global Land Evaporation Amsterdam Model and PML‐V2 datasets, but with some disagreement in the individual components, especially for evaporation. The surface soil moisture simulated by LPJ‐GUESS aligns with ESA‐CCI (v5.3) surface soil moisture datasets in most regions, with greatest discrepancies in subarctic areas. We attribute these discrepancies to two main sources: (1) absent or poor representation of processes such as river routing, storage and supply of water bodies, and cropland irrigation; and (2) uncertainties in both reference datasets and input to the model, including precipitation, soil texture, and land use.

Optimising CH4 simulations from the LPJ-GUESS model v4.1 using an adaptive Markov chain Monte Carlo algorithm

Abstract. The processes responsible for methane (CH4) emissions from boreal wetlands are complex; hence, their model representation is complicated by a large number of parameters and parameter uncertainties. The arctic-enabled dynamic global vegetation model LPJ-GUESS (Lund–Potsdam–Jena General Ecosystem Simulator) is one such model that allows quantification and understanding of the natural wetland CH4 fluxes at various scales, ranging from local to regional and global, but with several uncertainties. The model contains detailed descriptions of the CH4 production, oxidation, and transport controlled by several process parameters. Complexities in the underlying environmental processes, warming-driven alternative paths of meteorological phenomena, and changes in hydrological and vegetation conditions highlight the need for a calibrated and optimised version of LPJ-GUESS. In this study, we formulated the parameter calibration as a Bayesian problem, using knowledge of reasonable parameters values as priors. We then used an adaptive Metropolis–Hastings (MH)-based Markov chain Monte Carlo (MCMC) algorithm to improve predictions of CH4 emission by LPJ-GUESS and to quantify uncertainties. Application of this method on uncertain parameters allows for a greater search of their posterior distribution, leading to a more complete characterisation of the posterior distribution with a reduced risk of the sample impoverishment that can occur when using other optimisation methods. For assimilation, the analysis used flux measurement data gathered during the period from 2005 to 2014 from the Siikaneva wetlands in Southern Finland with an estimation of measurement uncertainties. The data are used to constrain the processes behind the CH4 dynamics, and the posterior covariance structures are used to explain how the parameters and the processes are related. To further support the conclusions, the CH4 flux and the other component fluxes associated with the flux are examined. The results demonstrate the robustness of MCMC methods to quantitatively assess the interrelationship between objective function choices, parameter identifiability, and data support. The experiment using real observations from Siikaneva resulted in a reduction in the root-mean-square error (RMSE), from 0.044 to 0.023 gC m−2 d−1, and a 93.89 % reduction in the cost function value. As a part of this work, knowledge about how CH4 data can constrain the parameters and processes is derived. Although the optimisation is performed based on a single site's flux data from Siikaneva, the algorithm is useful for larger-scale multi-site studies for a more robust calibration of LPJ-GUESS and similar models, and the results can highlight where model improvements are needed.

Tropical vegetation productivity and atmospheric methane over the last 40,000 years from model simulations and stalagmites in Sulawesi, Indonesia

AbstractRecent research has shown the potential of speleothem δ13C to record a range of environmental processes. Here, we report on 230Th-dated stalagmite δ13C records for southwest Sulawesi, Indonesia, over the last 40,000 yr to investigate the relationship between tropical vegetation productivity and atmospheric methane concentrations. We demonstrate that the Sulawesi stalagmite δ13C record is driven by changes in vegetation productivity and soil respiration and explore the link between soil respiration and tropical methane emissions using HadCM3 and the Sheffield Dynamic Global Vegetation Model. The model indicates that changes in soil respiration are primarily driven by changes in temperature and CO2, in line with our interpretation of stalagmite δ13C. In turn, modelled methane emissions are driven by soil respiration, providing a mechanism that links methane to stalagmite δ13C. This relationship is particularly strong during the last glaciation, indicating a key role for the tropics in controlling atmospheric methane when emissions from high-latitude boreal wetlands were suppressed. With further investigation, the link between δ13C in stalagmites and tropical methane could provide a low-latitude proxy complementary to polar ice core records to improve our understanding of the glacial–interglacial methane budget.

Mapping the global distribution of C4 vegetation using observations and optimality theory

Plants with the C4 photosynthesis pathway typically respond to climate change differently from more common C3-type plants, due to their distinct anatomical and biochemical characteristics. These different responses are expected to drive changes in global C4 and C3 vegetation distributions. However, current C4 vegetation distribution models may not predict this response as they do not capture multiple interacting factors and often lack observational constraints. Here, we used global observations of plant photosynthetic pathways, satellite remote sensing, and photosynthetic optimality theory to produce an observation-constrained global map of C4 vegetation. We find that global C4 vegetation coverage decreased from 17.7% to 17.1% of the land surface during 2001 to 2019. This was the net result of a reduction in C4 natural grass cover due to elevated CO2 favoring C3-type photosynthesis, and an increase in C4 crop cover, mainly from corn (maize) expansion. Using an emergent constraint approach, we estimated that C4 vegetation contributed 19.5% of global photosynthetic carbon assimilation, a value within the range of previous estimates (18–23%) but higher than the ensemble mean of dynamic global vegetation models (14 ± 13%; mean ± one standard deviation). Our study sheds insight on the critical and underappreciated role of C4 plants in the contemporary global carbon cycle.

Country-level estimates of gross and net carbon fluxes from land use, land-use change and forestry

Abstract. The reduction of CO2 emissions and the enhancement of CO2 removals related to land use are considered essential for future pathways towards net-zero emissions and mitigating climate change. With the growing pressure under global climate treaties, country-level land-use CO2 flux data are becoming increasingly important. So far, country-level estimates are mainly available through official country reports, such as the greenhouse gas inventories reported to the United Nations Framework Convention on Climate Change (UNFCCC). Recently, different modelling approaches, namely dynamic global vegetation models (DGVMs) and bookkeeping models, have moved to higher spatial resolutions, which makes it possible to obtain model-based country-level estimates that are globally consistent in their methodology. To progress towards a largely independent assessment of country reports using models, we analyse the robustness of country-level CO2 flux estimates from different modelling approaches in the period 1950–2021 and compare them with estimates from country reports. Our results highlight the general ability of modelling approaches to estimate land-use CO2 fluxes at the country level and at higher spatial resolution. Modelled land-use CO2 flux estimates generally agree well, but the investigation of multiple DGVMs and bookkeeping models reveals that the robustness of their estimates strongly varies across countries, and substantial uncertainties remain, even for top emitters. Similarly, modelled land-use CO2 flux estimates and country-report-based estimates agree reasonably well in many countries once their differing definitions are accounted for, although differences remain in some other countries. A separate analysis of CO2 emissions and removals from land use using bookkeeping models also shows that historical peaks in net fluxes stem from emission peaks in most countries, whereas the long-term trends are more connected to removal dynamics. The ratio of the net flux to the sum of CO2 emissions and removals from land use (the net-to-gross flux ratio) underlines the spatio-temporal heterogeneity in the drivers of net land-use CO2 flux trends. In many tropical regions, net-to-gross flux ratios of about 50 % are due to much larger emissions than removals; in many temperate countries, ratios close to zero show that emissions and removals largely offset each other. Considering only the net flux thus potentially masks large emissions and removals and the different timescales upon which they act, particularly if averaged over countries or larger regions, highlighting the need for future studies to focus more on the gross fluxes. Data from this study are openly available via the Zenodo portal at https://doi.org/10.5281/zenodo.8144174 (Obermeier et al., 2023).

Connecting competitor, stress-tolerator and ruderal (CSR) theory and Lund Potsdam Jena managed Land 5 (LPJmL 5) to assess the role of environmental conditions, management and functional diversity for grassland ecosystem functions

Abstract. Forage offtake, leaf biomass and soil organic carbon storage are important ecosystem services of permanent grasslands, which are determined by climatic conditions, management and functional diversity. However, functional diversity is not independent of climate and management, and it is important to understand the role of functional diversity and these dependencies for ecosystem services of permanent grasslands, since functional diversity may play a key role in mediating impacts of changing conditions. Large-scale ecosystem models are used to assess ecosystem functions within a consistent framework for multiple climate and management scenarios. However, large-scale models of permanent grasslands rarely consider functional diversity. We implemented a representation of functional diversity based on the competitor, stress-tolerator and ruderal (CSR) theory and the global spectrum of plant form and function into the Lund Potsdam Jena managed Land (LPJmL) dynamic global vegetation model (DGVM) forming LPJmL-CSR. Using a Bayesian calibration method, we parameterised new plant functional types (PFTs) and used these to assess forage offtake, leaf biomass, soil organic carbon storage and community composition of three permanent grassland sites. These are a temperate grassland and a hot and a cold steppe for which we simulated several management scenarios with different defoliation intensities and resource limitations. LPJmL-CSR captured the grassland dynamics well under observed conditions and showed improved results for forage offtake, leaf biomass and/or soil organic carbon (SOC) compared to the original LPJmL 5 version at the three grassland sites. Furthermore, LPJmL-CSR was able to reproduce the trade-offs associated with the global spectrum of plant form and function, and similar strategies emerged independent of the site-specific conditions (e.g. the C and R PFTs were more resource exploitative than the S PFT). Under different resource limitations, we observed a shift in the community composition. At the hot steppe, for example, irrigation led to a more balanced community composition with similar C, S and R PFT shares of aboveground biomass. Our results show that LPJmL-CSR allows for explicit analysis of the adaptation of grassland vegetation to changing conditions while explicitly considering functional diversity. The implemented mechanisms and trade-offs are universally applicable, paving the way for large-scale application. Applying LPJmL-CSR for different climate change and functional diversity scenarios may generate a range of future grassland productivities.

Atmospheric water demand constrains net ecosystem production in subtropical mangrove forests

Subtropical mangroves have great potential to sequester atmospheric carbon dioxide (CO2) and thus significantly contribute to climate change mitigation. Meanwhile, the carbon cycling processes in subtropical mangroves are vulnerable to the changing climatic conditions, such as the warming-induced vapor pressure deficit (VPD) increases. However, the impacts of VPD on net ecosystem production (NEP) in subtropical mangroves remain poorly understood due to a lack of detailed assessment of NEP responses to VPD among different mangrove forests over a long-term observational period. In this study, we deployed eddy-covariance systems to measure net ecosystem CO2 exchange at three subtropical mangrove forests for 16 site-years, comprising two natural and one artificial forests. We employed a state-of-the-art data-driven modeling approach (i.e., Shapley additive explanations (SHAP) framework based on extreme gradient boosting model), which enabled us to explore the interactive effects of meteorological and tidal factors (e.g. salinity) with VPD on the mangrove NEP. We revealed that air temperature, global solar radiation (GR) and wind speed have significant interactions on the response of NEP to VPD stress in subtropical mangroves. For instance, when GR was high, the SHAP interaction values of VPD and GR on NEP decreased with increasing VPD, but when GR was low, the trends were the opposite. However, instead of identifying interactive effects between tidal salinity and VPD on mangrove NEP, we came across potential independent influence of salinity on the same. SHAP analysis was also able to disentangle the impact of VPD from other abiotic drivers. Thus, we evaluated the threshold effect of VPD stress on NEP loss in subtropical mangroves and observed a range of 2.50–2.95 kPa. Above this range, VPD stress leveled off. The subtropical mangrove responses to VPD should be therefore considered in the dynamic global vegetation models to increase the accuracy in carbon cycle simulations in the future.

African rainforest moisture contribution to continental agricultural water consumption

Precipitation is essential for food production in Sub-Saharan Africa, where more than 80 % of agriculture is rainfed. Although ∼40 % of precipitation in certain regions is recycled moisture from Africa's tropical rainforest, there needs to be more knowledge about how this moisture supports the continent's agriculture. In this study, we quantify all moisture sources for agrarian precipitation (African agricultural precipitationshed), the estimates of African rainforest's moisture contribution to agricultural precipitation, and the evaporation from agricultural land across the continent. Applying a moisture tracking model (UTRACK) and a dynamic global vegetation model (LPJmL), we find that the Congo rainforest (>60 % tree cover) is a crucial moisture source for many agricultural regions. Although most of the rainforest acreage is in the DRC, many neighboring nations rely significantly on rainforest moisture for their rainfed agriculture, and even in remote places, rainforest moisture accounts for ∼10–20 % of agricultural water use. Given continuous deforestation and climate change, which impact rainforest areas and resilience, more robust governance for conserving the Congo rainforest is necessary to ensure future food production across multiple Sub-Saharan African countries.

Water supply following drought: Effects on drougt legacy and resilience in a tropical forest - A case study in Xishuangbanna, China

Climate change has made it essential to understand the effects of drought on carbon dynamics in tropical forest ecosystems. However, research on the postdrought (PD) period is relatively limited, particularly with regard to the resilience of the carbon cycle affected by available water. To determine the effects of water availability during PD on resilience and legacy effects in tropical forests, we used a spatially explicit individual-based dynamic global vegetation model (SEIB-DGVM) to simulate gross primary productivity (GPP) to proxy the process. We established multiple water supply scenarios based on ecological water demands and regional climatic standard precipitation through the composite of gradients in quantity and duration to characterize variations in environmental water availability after droughts end. Under different PD water supply patterns, the degree of GPP resilience and its legacy effects varied significantly. Specifically, ecosystem resilience to drought was positively correlated with the percentage of PD water supply and negatively correlated with the duration until water conditions returned to normal, while the negative legacy effect of drought on GPP exhibited the opposite relationship. After a drought, a rapid high-level water supply promoted stronger GPP resiliency and moderated the drought legacy effects. Although legacy effects varied with water supply patterns, the negative legacy effect on GPP appeared only in the first year after drought. Observing biomass changes in roots, leaves, woody plants, and herbaceous plants after droughts, it is clear that tropical forest ecosystems can adjust their carbon allocation strategies in response to differences in water availability. To more accurately predict the fate of tropical forests in response to climate change, it is necessary to assess the impact of PD environmental climate condition changes on their dynamic resilience processes.

Temporal variability of observed and simulated gross primary productivity, modulated by vegetation state and hydrometeorological drivers

Abstract. The gross primary production (GPP) of the terrestrial biosphere is a key source of variability in the global carbon cycle. It is modulated by hydrometeorological drivers (i.e. short-wave radiation, air temperature, vapour pressure deficit and soil moisture) and the vegetation state (i.e. canopy greenness, leaf area index) at instantaneous to interannual timescales. In this study, we set out to evaluate the ability of GPP models to capture this variability. Eleven models were considered, which rely purely on remote sensing data (RS-driven), meteorological data (meteo-driven, e.g. dynamic global vegetation models; DGVMs) or a combination of both (hybrid, e.g. light-use efficiency, LUE, models). They were evaluated using in situ observations at 61 eddy covariance sites, covering a broad range of herbaceous and forest biomes. The results illustrated how the determinant of temporal variability shifts from meteorological variables at sub-seasonal timescales to biophysical variables at seasonal and interannual timescales. RS-driven models lacked the sensitivity to the dominant drivers at short timescales (i.e. short-wave radiation and vapour pressure deficit) and failed to capture the decoupling of photosynthesis and canopy greenness (e.g. in evergreen forests). Conversely, meteo-driven models accurately captured the variability across timescales, despite the challenges in the prognostic simulation of the vegetation state. The largest errors were found in water-limited sites, where the accuracy of the soil moisture dynamics determines the quality of the GPP estimates. In arid herbaceous sites, canopy greenness and photosynthesis were more tightly coupled, resulting in improved results with RS-driven models. Hybrid models capitalized on the combination of RS observations and meteorological information. LUE models were among the most accurate models to monitor GPP across all biomes, despite their simple architecture. Overall, we conclude that the combination of meteorological drivers and remote sensing observations is required to yield an accurate reproduction of the spatio-temporal variability of GPP. To further advance the performance of DGVMs, improvements in the soil moisture dynamics and vegetation evolution are needed.

Global Carbon Budget 2023

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO2 products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2022, EFOS increased by 0.9 % relative to 2021, with fossil emissions at 9.9±0.5 Gt C yr−1 (10.2±0.5 Gt C yr−1 when the cement carbonation sink is not included), and ELUC was 1.2±0.7 Gt C yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1±0.8 Gt C yr−1 (40.7±3.2 Gt CO2 yr−1). Also, for 2022, GATM was 4.6±0.2 Gt C yr−1 (2.18±0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.8±0.4 Gt C yr−1, and SLAND was 3.8±0.8 Gt C yr−1, with a BIM of −0.1 Gt C yr−1 (i.e. total estimated sources marginally too low or sinks marginally too high). The global atmospheric CO2 concentration averaged over 2022 reached 417.1±0.1 ppm. Preliminary data for 2023 suggest an increase in EFOS relative to 2022 of +1.1 % (0.0 % to 2.1 %) globally and atmospheric CO2 concentration reaching 419.3 ppm, 51 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2022, with a near-zero overall budget imbalance, although discrepancies of up to around 1 Gt C yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living-data update documents changes in methods and data sets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2023 (Friedlingstein et al., 2023).

Climate and dispersal limitation drive tree species range shifts in post-glacial Europe: results from dynamic simulations

IntroductionThe ability of species to colonize newly suitable habitats under rapid climate change can be constrained by migration processes, resulting in a shift of the leading edge lagging behind the ameliorating climate, i.e. migration lag. The importance and extent of such migration lags during the forest expansion after the Last Glacial Maximum (LGM) are still debated. Similarly, the relative importance of the main drivers of post-glacial vegetation dynamics (temperature, dispersal limitation, and competition) is still discussed in the literature.MethodsWe used the dynamic global vegetation model LPJ-GM 2.0 to reconstruct the range shifts of 16 competing major European tree species after the LGM (18.5 ka BP) until recent times (0 ka BP). We simulated two dispersal modes by allowing free establishment whenever the climatic conditions suited the species (free dispersal), or by accounting for migration processes in the simulated vegetation dynamics (dispersal limitation). We then calculated thermal and range shift velocities, competition at establishment, thermal and dispersal lags for each species and dispersal mode. Finally, we compared our simulated range shift velocities with pollen-derived migration rates.ResultsThe simulation assuming limited dispersal resulted in more accurate migration rates as compared to pollen-derived migration rates and spreading patterns. We found no marked migration lags in the post-glacial establishment of pioneer species (Pinus sylvestris and Betula pubescens). Under the free dispersal mode, the remaining temperate species expanded rapidly and almost synchronously across central Europe upon climate warming (Bølling-Allerød interstadial). Differently, the northward spread of temperate species simulated under dispersal limitation happened mainly during the Holocene and in successive waves, with late spreaders (e.g. Fraxinus excelsior) experiencing multi-millennial dispersal lags and higher competition.DiscussionOur simulation under dispersal constraints suggests that the post-glacial tree expansion in Europe was mainly driven by species-specific thermal requirements and dispersal capacity, which in turn affected the order of taxa establishment and thus the degree of competition. Namely, taxa with less cold-tolerance and relatively low dispersal ability experienced the highest migration lags, whereas the establishment of pioneer species was mostly in equilibrium with the climate.

Biome classification influences current and projected future biome distributions

AbstractAimBiome classification schemes are widely used to map biogeographic patterns of vegetation formations on large spatial scales. Future climate change will influence biome patterns, and vegetation models can be used to assess the susceptibility of biomes to experience transitions. However, biome classification is not unique, and various classification schemes and biome maps exist. Here, we aimed to assess how the choice of biome classification schemes influences current and projected future biome patterns.LocationAfrica, Australia, Tropical Asia.Time period2000–2099.Major taxa studiedTropical vegetation.MethodsWe used the adaptive dynamic global vegetation model version 2 (aDGVM2) to simulate vegetation in the study region. We classified vegetation into biomes using (1) a classification scheme based on the cover of functional types, (2) a cluster analysis based on the cover of functional types and (3) a cluster analysis based on trait patterns simulated by the aDGVM2. We compared the resulting biome maps to multiple observation‐based biome maps and quantified differences in projected biome changes under the RCP8.5 scenario for the different classification schemes.ResultsAs expected, biome patterns were strongly related to the scheme used for biome classification. The highest data‐model agreement was derived for a cluster analysis using 21 simulated traits. Traits related to size were most important for classification. Considering all classification schemes, the area projected to undergo biome transitions under climate change varied between 16.5% and 32.1%. Despite this variability, different schemes consistently showed that grassland and savanna areas are most susceptible to climate change, whereas tropical forests and deserts are stable. Our results demonstrate that traits simulated by aDGVM2 are appropriate to delimit biomes.Main conclusionsStudies projecting biome patterns and transitions under current and future climate should consider applying different biome classification schemes to avoid biases in such projections caused by biome classification schemes.

Trend in seasonal amplitude of northern net ecosystem production: Simulated results from IAP DGVM in CAS-ESM2

Net ecosystem production (NEP) over northern latitudes has an obvious seasonal variation, which dominantly drives the seasonal variation in atmospheric CO2. Investigating these carbon processes and understanding the underlying drivers is a key issue of climate research. This study focuses on the seasonal amplitude of northern NEP (40°–90°N) and investigates its trend in the period 1990–2014 by using the dynamic global vegetation model (IAP DGVM) in the second version of the Chinese Academy of Sciences Earth System Model (CAS-ESM2). Basing on a spin-up simulation, the authors conduct a control simulation to evaluate the modeled trend in the seasonal amplitude of northern NEP, and three sensitivity simulations to detect the contributions of the climate and atmospheric CO2. The results show that the modeled seasonal amplitude of northern NEP increases significantly from 1990 to 2014, with a trend of 9.69 TgC month−1 yr−1, which is mainly because of the increasing maximum NEP. The positive trend is largely reduced when CO2 fertilization effects and climatic effects are separately excluded. These significant reductions suggest important effects of atmospheric CO2 and climate change on the seasonal amplitude of northern NEP. Even though model uncertainties remain, the results favor further development of IAP DGVM in accurately simulating the terrestrial carbon cycle, and also provide an important reference for the application of CAS-ESM in further investigating carbon–climate interactions.摘要北方陆地净生态系统生产力(NEP)具有明显的季节变化特征, 这是大气CO2季节变化的关键驱动. 研究这些碳循环过程并理解潜在的驱动因素是气候研究的一个关键问题. 本文利用第二代中国科学院地球系统模式(CAS-ESM2)中的全球植被动态模型(IAP DGVM), 研究了1990−2014年北方NEP(40°−90°N)的季节振幅及其变化趋势. 在初始化试验的基础上, 本文开展了一个控制试验来评估模拟的北方NEP季节幅度的变化趋势, 同时开展了三个敏感性试验来研究气候和大气CO2的贡献. 结果表明: 1990−2014年, 模拟的北方NEP季节振幅显著增加, 趋势为9.69 万吨碳/月/年, 这主要是由于最大NEP增加所致. 当分别排除CO2施肥效应和气候效应时, 上述增加趋势大大减弱.这些显著的减少表明大气CO2和气候变化对北方NEP的季节性振幅有重要影响. 尽管模式存在不确定性, 但这些结果有利于进一步提升IAP DGVM对陆地碳循环的精确模拟, 也为CAS-ESM研究碳-气候相互作用的应用提供了重要参考.

A theory of demographic optimality in forests

Carbon uptake by the land is a key determinant of future climate change. Unfortunately, Dynamic Global Vegetation Models have many unknown internal parameters which leads to significant uncertainty in projections of the future land carbon sink. By contrast, observed forest inventories in both Amazonia and the USA show strikingly common tree-size distributions, pointing to a simpler modelling paradigm. The curvature of these size-distributions is related to the ratio of mortality to growth in Demographic Equilibrium Theory (DET). We extend DET to include recruitment limited by competitive exclusion from existing trees. From this, we find simultaneous maxima of tree density and biomass in terms of respectively the ratio of mortality to growth and the proportion of primary productivity allocated to reproduction, an idea we call Demographic Optimality (DO). Combining DO with the ratio of mortality to growth common to the US and Amazon forests, results in the prediction that about an eighth of productivity should be allocated to reproduction, which is broadly consistent with observations. Another prediction of the model is that seed mortality should decrease with increasing seed size, such that the advantage of having many small seeds is nullified by the higher seed mortality. Demographic Optimality is therefore consistent with the common shape of tree-size distributions seen in very different forests, and an allocation to reproduction that is independent of seed size.